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Abstract 

A search for single-top production, ep —>■ etX, has been made with the ZEUS 
detector at HERA using an integrated luminosity of 130.1 pb~^. Events from 
both the leptonic and hadronic decay channels of the W boson resulting from the 
decay of the top quark were sought. For the leptonic mode, the search was made 
for events with isolated high-energy leptons and significant missing transverse 
momentum. For the hadronic decay mode, three-jet events in which two of the 
jets had an invariant mass consistent with that of the W were selected. No 
evidence for top production was found. The results are used to constrain single- 
top production via flavour-changing neutral current (FCNC) transitions. The 
ZEUS limit excludes a substantial region in the FCNC twy coupling not ruled 
out by other experiments. 
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1 Introduction 



The Standard Model (SM) of the fundamental interactions presently provides an accurate 
description of the phenomena observed in both low- and high-energy reactions of elemen- 
tary particles. As probes in search for physics beyond the SM, observables sensitive to 
flavour-changing neutral current (FCNC) interactions are particularly useful, since the 
SM rates are very small due to the GIM mechanism [1]. The FCNC interactions involving 
the top quark [2,3], which has a mass of the order of the electroweak energy scale, offer 
a potentially new view of physics beyond the SM. 

The FCNC-induced couplings of the type tuV or tcV (with y = 7, Z°) have been explored 
in pp collisions at the Tevatron by searching for the top-quark decays t uV and t — 
cV [4]. The same couplings involving the top quark were investigated in e^e~ interactions 
at LEP2 by searching for single-top production through the reactions e+e~ — tu (+c.c.) 
and e+e~ — > tc (+c.c.) [5,6]. No evidence for such interactions was found at either 
accelerator and limits were set on the branching ratios B{t 57) and B{t ^ qZ). 

In ep collisions at the HERA collider, top quarks can only be singly produced. In the SM, 
single-top production proceeds through the charged current (CC) reaction ep — > utbX [7]. 
Since the SM cross section at HERA is less than 1 fb [8], any observed single-top event 
in the present data can be attributed to physics beyond the SM. The FCNC couplings, 
tuV or tcV, would induce the neutral current (NC) reaction ep etX [3,9], in which 
the incoming lepton exchanges a 7 or Z with an up-type quark in the proton, yielding 
a top quark in the final state. Due to the large Z mass, this process is most sensitive 
to a coupling of the type tq'y. Furthermore, large values of x, the fraction of the proton 
momentum carried by the struck quark, are needed to produce a top. Since the w-quark 
parton distribution function (PDF) of the proton is dominant at large x, the production 
of single top quarks is most sensitive to a coupling of the type twj (see Fig. [T]). 



2 Theoretical framework 

Deviations from the SM predictions due to FCNC transitions involving the top quark can 
be parameterised in terms of couplings of the type tuV (with V = Z°) and described 
by an effective Lagrangian of the form [10] 

A^eflF = e ett — titu-y u A^" + — t 7^ vtuz u + h.c, (1) 

where e (et) is the electron (top-quark) electric charge, g is the weak coupling constant, 
6w is the weak mixing angle, 0"^,^ = |(7'^7'^ — 7'^7^), A is an effective cutoff which, by 
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convention, is set to the mass of the top quark, M^op, taken as 175 GeV, q is the momentum 
of the gauge boson and {Z'^) is the photon {Z) field. In the following, it was assumed 
that the magnetic coupling and the vector coupling Vtuz are real and positive. The 
values of Ktu-y and Vtuz in the SM are zero at tree level and extremely small at the one-loop 
level. 

The cross section for the process ep —>■ etX was calculated as a function of Ktu-y includ- 
ing next-to- leading-order (NLO) QCD corrections in the eikonal approximation [9]. The 
renormalisation (/ir) and factorisation (fip) scales were chosen to he = fip = ^top- 
The strong coupling constant, Og, was calculated at two loops with Ai^ = 220 MeV, cor- 
responding to as{Mz) = 0.1175. The calculations were performed using the MRST99 [11] 
parameterisations of the proton PDFs. The uncertainty of the results due to terms be- 
yond NLO, estimated by varying = between Mtop/2 and 2Mtop, was tl'^Yo (^1^%) 
at a centre-of-mass energy of 318 (300) GeV. The uncertainties of the results due to that 
on as{Mz) and on the proton PDFs were ±2% and ±4%, respectively. The variation of 
the cross section on Mtop was approximately ±20% (±25%) for AMtop = ±5 GeV at a 
centre-of-mass energy of 318 (300) GeV. 



3 Experimental conditions 

The data samples were collected with the ZEUS detector at HERA and correspond to 
an integrated luminosity of 47.9 ± 0.9 (65.5 ± 1.5) pb~^ for e'^p collisions taken during 
1994-1997 (1999-2000) and 16.7±0.3 pb'^ for e'p colhsions taken during 1998-99. During 
1994-97 (1998-2000), HERA operated with protons of energy Ep = 820 GeV (920 GeV) 
and positrons or electrons of energy = 27.5 GeV, yielding a centre-of-mass energy of 
^ = 300 GeV (318 GeV). 

The ZEUS detector is described in detail elsewhere [12,13]. The main components used in 
the present analysis were the central tracking detector (CTD) [14], positioned in a 1.43 T 
solenoidal magnetic field, and the uranium-scintillator sampling calorimeter (CAL) [15]. 

Tracking information is provided by the CTD, in which the momenta of tracks in the 
polar-angle^ region 15° < 6 < 164° are reconstructed. The CTD consists of 72 cylindrical 
drift chamber layers, organised in nine superlayers. The relative transverse momentum, 
Pt, resolution for full-length tracks can be parameterised as <j{pt)/pt = 0.0058 pt © 
0.0065 © 0.0014/pr, with pr in GeV. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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The CAL covers 99.7% of the total sohd angle. It is divided into three parts with a 
corresponding division in 6', as viewed from the nominal interaction point: forward (FCAL, 
2.6° < e < 36.7°), barrel (BCAL, 36.7° <e < 129.1°), and rear (RCAL, 129.1° < 6 < 
176.2°). Each of the CAL parts is subdivided into towers which in turn are segmented 
longitudinally into one electromagnetic (EMC) and one (RCAL) or two (FCAL, BCAL) 
hadronic (HAC) sections. The smallest subdivision of the CAL is called a cell. Under 
test-beam conditions, the CAL single-particle energy resolution is a{E)/E = 18%/ 
for electrons and a{E)/E = 35% /VE for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — ^ e'jp. The resulting 
small-angle energetic photons were measured by the luminosity monitor [16], a lead- 
scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3.1 Trigger conditions 

A three-level trigger was used to select events online [12, 17]. At the first level, events 
were selected using criteria based on either the transverse energy or missing transverse 
momentum measured in the CAL. Events were accepted with a low threshold on these 
quantities when a coincidence with CTD tracks from the event vertex was required, while 
a higher threshold was used for events with no CTD tracks. 

At the second level, timing information from the CAL was used to reject events incon- 
sistent with an ep interaction. In addition, the topology of the CAL energy deposits 
was used to reject non-ep background events. Cuts on the missing transverse momentum 
of 6 GeV (9 GeV for events without CTD tracks) or on the total transverse energy of 
8 GeV, excluding the eight CAL towers immediately surrounding the forward beampipe, 
were applied. 

At the third level, track reconstruction and vertex finding were performed and used to 
reject events with a vertex inconsistent with the distribution of ep interactions. Events 
with missing transverse momentum in excess of 7 GeV or containing at least two jets 
with transverse energy i?;^* > 6 GeV and pseudorapidity t]^^^ < 2.5 were accepted; the 
latter condition was based upon the application of a jet-finding cone algorithm with radius 
-R = 1 applied to the CAL cell energies and positions. 

4 Monte Carlo simulation 

Samples of events were generated using Monte Carlo (MC) simulations to determine the 
selection efficiency for the signal of single-top production through FCNC processes and to 
estimate background rates from SM processes. The generated events were passed through 



3 



the GEANT 3.13-based [18] ZEUS detector- and trigger-simulation programs [12]. They 
were reconstructed and analysed by the same program chain as the data. 

Single-top production through FCNC processes in ep collisions was simulated using the 
HEXF generator [19]. Samples of events were generated assuming top-quark masses of 
170, 175 and 180 GeV. Initial-state radiation from the lepton beam was included using 
the Weizsacker- Williams approximation [20] . The hadronic final state was simulated using 
the matrix-element and parton-shower model of LEPTO [21] for the QCD cascade and 
the Lund string model [22] as implemented in JETSET [23] for the hadronisation. The 
MRSA [24] parameterisations of the proton PDFs were used. 

The most important background to the positron-decay channel of the W in the chain 
t — *■ bW'^ — > be'^u arose from NC deep inelastic scattering (DIS). Two-photon processes 
provide a source of high-p-r leptons that were a significant background to the muon-decay 
channel of the W in the chain t hW^ bfi'^u. In addition, single- production was 
a significant source of background to t ^ bW~^, in both the positron- and muon-decay 
channels of the W. The dominant source of background for the hadronic-decay channel 
of W in the chain t — > bW~^ — > bqq' was multi-jet production from QCD processes. 

Several MC programs were used to simulate the different background processes. The 
NC DIS events were generated using the LEPTO 6.5 program [21] interfaced to HER- 
ACLES 4.6.1 [25] via DJANGOH 1.1 [26]. The HERACLES program includes photon 
and Z exchanges and first-order electroweak radiative corrections. The QCD cascade was 
modelled with the colour-dipole model [27] by using the ARIADNE 4.08 program [28] 
and including the boson-gluon-fusion process. As an alternative, samples of events were 
generated using the model of LEPTO based on first-order QCD matrix elements plus 
parton showers (MEPS). In both cases, the hadronisation was performed using the Lund 
string model. The CTEQ5D [29] parameterisations for the proton PDFs were used. Two- 
photon processes were simulated using the generator GRAPE 1.1 [30], which includes 
dilepton production via 77, Z7 and ZZ processes and considers both elastic and inelastic 
production at the proton vertex. Single-l^ production was simulated using the event 
generator EPVEC [31], which did not include hard QCD radiation. Recent cross-section 
calculations including higher-order QCD corrections [32] and using the CTEQ4M [33] 
(ACFGP [34]) proton (photon) PDFs were used to reweight the EPVEC event samples. 
Multi-jet QCD production at low Q^, where is the virtuality of the exchanged photon, 
was simulated using PYTHIA 5.7 [35]. In this generator, the partonic processes were 
simulated using leading-order (LO) matrix elements, with the inclusion of initial- and 
final-state parton showers. Hadronisation was performed using the Lund string model. 
The MRSA (GRV-HO [36]) parameterisations of the proton (photon) PDFs were used. 
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5 Signatures of FCNC-induced single-top production 



Single-top production via the FCNC coupling at the twy vertex in ep collisions at HERA, 
ep — > etX, is predicted to proceed predominantly through the exchange of a quasi- real 
photon between the beam electron or positron and a valence u quark in the proton (see 
Fig. P). According to the signal MC simulation, the scattered electron or positron escapes 
through the rear beampipe, outside the CAL acceptance, in 65% of the events. 

In this analysis, the top-quark search was optimised for the decay t bW~^. In the 
leptonic decay channel of the W, the signal for such events is the presence of an isolated 
high-energy lepton, significant missing transverse momentum arising from the emitted 
neutrino and a jet stemming from the 6-quark decay. In the hadronic decay channel of 
the W, the signal is the presence of three jets in the final state with the dijet invariant- 
mass distribution for the correct pair of jets peaking at the mass of the W boson, Mw, 
and the three-jet invariant-mass distribution peaking at Mtop. 



6 Leptonic channel 
6.1 Data selection 

Events with isolated high-energy leptons (e^ or fi"^), significant missing transverse mo- 
mentum and a jet were selected. Similar previous analyses have been done by the HI [37] 
and ZEUS [38] Collaborations. Positron candidates were identified using an algorithm 
that combined CAL and CTD information [39] . Muons were identified by the coincidence 
of a track in the CTD with significant transverse momentum and CAL energy deposits 
consistent with those expected from a minimum ionizing particle. The charge informa- 
tion on the candidates was not used and they are generically referred to as positrons and 
muons. The main selection criteria are: 

• cuts on the CAL timing and Z coordinate {\Z\ < 50 cm) of the event vertex and 
algorithms based on the pattern of tracks in the CTD were used to reject events not 
originating from ep collisions; 

• the track associated with the positron or muon candidate was required to have p^'^'^^ > 
5 GeV. To reduce the NC DIS background, the track was required to have 6 < 115°. 
In addition, it must have passed through at least three radial superlayers of the CTD 
(corresponding to 9 > 17°) and be isolated. Two isolation variables were defined for 
a given track using the separation R in the rj — (p plane, where R = (Arj)'^ + {A(py. 
The variable Dj^t was defined as the distance from the nearest jet axis, while -Dtrack 
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was the distance from the nearest neighbouring track in the event. Events containing 
tracks with Djct > 1 and -Dtrack > 0.5 were selected; 

• p^^^ > 20 GeV, where p^^^ is the missing transverse momentum as measured with the 
CAL. It was reconstructed using the energy deposited in the CAL cells, after correc- 
tions for non- uniformity and dead material located in front of the CAL [40]. Energy 
deposits originating from identified muons were excluded from the measurement of 

„CAL. 

Pt ) 

• the presence of at least one jet with transverse energy i?^* above 5 GeV and — 1 < 
77-^'^* < 2.5 was required. The longitudinally invariant kx cluster algorithm [41] was 
used in the inclusive mode [42] to reconstruct jets from the energy deposits in the 
CAL cells. The jet search was performed in the f] — plane of the laboratory frame. 
The axis of each jet was defined according to the Snowmass convention [43], where 
j^jet (^(^jct^ ^g^g ^]-^g transverse-energy-weighted mean pseudorapidity (azimuth) of all 
the cells belonging to that jet. The jet transverse energy was reconstructed as the 
sum of the transverse energies of the cells belonging to the jet and was corrected for 
detector effects such as energy losses in the inactive material in front of the CAL [44]. 
In the leptonic channel, only those jets for which the electromagnetic-energy fraction 
was below 0.9 and Rqo% > 0.1, where Rgo% is the radius of the cone in the r] — ip plane 
concentric to the jet axis that contains 90% of the jet energy, were considered; 

• in events with an identified positron candidate, the acoplanarity angle, $acop5 was 
defined as the azimuthal separation of the outgoing positron and the vector in the 
{X, y)-plane that balances the hadronic system. For well measured NC DIS events, 
the acoplanarity angle is close to zero, while a large $acop indicates large missing 
energy, as expected from top-quark decays. To reduce the background from NC DIS 
processes, the acoplanarity angle was required to be greater than 8°. 

The selected data sample contained 36 events, 24 of which had a positron candidate and 
12 a muon candidate (see Table H}. 



6.2 Comparison with Monte Carlo simulations 

The properties of the selected events were studied in detail and compared with the MC 
predictions of the SM. Figures |2ti')-c) show the acoplanarity, the transverse momentum 
of the hadronic system, p^'^, and the transverse momentum of the positron candidate 
as measured in the CAL, p^, for those events with an identified positron candidate. 
Figures |21i)-f ) show the CAL transverse momentum corrected for the muon momentum 
measured by the CTD, p^^ = \/ {Px^^ + Px^ + {Py^^ + PyVy Pr^ transverse 
momentum of the track associated with the muon candidate, p^, for events with an 
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identified muon candidate. In each case, the distribution of data events can be accounted 
for by the simulation of SM processes. The SM expectation for the positron channel is 
dominated by NC DIS and that for the muon channel by two-photon processes. 

6.3 Results of the search in the leptonic channel 

The final event selection for the blu final state required a high-p^ jet and missing en- 
ergy. The cuts were optimised using the simulations of both the SM background and the 
expected single-top signal. The cuts used were: 

• p^'^ > 40 GeV for both the positron and muon decays; 

• 6 = cos 9i) = ~Pz)i < 47 GeV for the positron decay, where the sum 
runs over all GAL energy deposits with corrected energy Ei and polar angle 6i [40]. 
For fully contained NG DIS events, 6 peaks at 55 GeV, i.e. twice the lepton beam 
energy, which follows from energy-momentum conservation; 

• > 10 GeV for the muon decay. 

After applying these requirements, no event remained in the data sample. The efficiency 
for detecting single-top production in the leptonic channel was 34% for the positron decay 
and 33% for the muon decay. These efficiencies do not include the branching ratio of the 
top-quark decay in the corresponding channel. 

In a recent study [45], the HI GoUaboration has reported an excess of events for p^'^ > 
25 GeV. The number of selected events in each channel with p^'^ > 25 GeV for the present 
analysis is also listed in Table ^ These results are in agreement with the expectations 
from the SM. 

7 Hadronic channel 

The data used for this channel correspond to a slightly reduced luminosity of 127.2 pb~^. 
7.1 Data selection 

The expected signature for the hadronic-decay channel of single-top production through 
the FGNG twy coupling is three jets with large E^r^^ and no significant missing transverse 
momentum. Since it is expected that for the bulk of the events the scattered positron 
escapes through the rear beam pipe, NG DIS events with > 1 GeV^ were rejected. The 
data selection used similar criteria as reported in a previous publication [44]. Jets were 



7 



found in the hadronic final state using the same algorithm as described in Section 16.11 
The main selection criteria are: 

• cuts on the Z coordinate (—38 < Z < 32 cm) of the event vertex, the number of 
tracks pointing to the vertex and the number of tracks compatible with an interaction 
upstream in the direction of the proton beam were used to reject events not originating 
from ep collisions; 

• the presence of at least three jets within the pseudorapidity range —1 < 77-''^* < 2.5 
was required. The three highest-E'^?* jets in the event, ordered according to decreasing 

were further required to satisfy > 40, 25, 14 GeV; 

• CC DIS events were rejected by requiring the missing transverse momentum to be 
small compared to the total transverse energy, E*?*, i.e. p'^^^ / ^/E^ < 2 VGeV; 

• NC DIS events with an identified scattered-positron candidate [46] in the CAL were 
removed from the sample using the method described in an earlier publication [47]; 

• 8.8 < 6 < 52.2 GeV. The upper cut removed unidentified NC DIS events and the lower 
cut rejected proton beam-gas interactions. 

The selected sample contained 348 events. 

The invariant mass of jets k and / was determined using the corrected jet transverse 
energies, as explained in Section lUTTl and jet angular variables according to the formula 

ikfjj = ^2E}^^'^E}^^'\cosh (r/J^t.^ - rji^^^) - cos {<^i^''^ - (^j^*'')]. 
The three-jet invariant mass, was reconstructed using the formula 

M^j = 2£;j.^*'''Ej?*'^[cosh (r/j^t'k _ ^jct,i) _ ^os ((^J^t'^ _ ^^jct,!)]^ 

y k<i 

where the sum runs over k,l = 1,2,3. The average resolution in M^^ was 8% for M^j > 
50 GeV and the distribution of M-^j for all pairs of jets in a sample of MC signal events is 
shown in Fig. ^.). The average resolution in M^j was 4% for M^j > 80 GeV and the M^j 
distribution in a sample of MC signal events is shown in Fig. IHb)- Cuts on M-^j and M^j 
were used to search for a signal of single-top production in the hadronic channel. 

7.2 Comparison with Monte Carlo simulations 

The properties of the selected events were studied in detail and were compared with 
the MC predictions of the SM processes. The MC distributions were normalised to the 
number of events in the data with M^^ < 159 GeV, i.e. outside the region where the signal 
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for single-top production is expected. The resulting normalisation factor was 1.11 ±0.08, 
which can be attributed to higher-order QCD corrections to the jet cross sections. The 
simulations of SM processes provide a reasonable description of the E^r^^ and rj^^^ data 
distributions (not shown). The distribution of M^^ for all pairs of jets in an event and 
that of M^j are presented in Figs. Et) and b), respectively. The M^j distribution shows a 
steep fall-off from M^j ~ 130 GeV to 240 GeV. The SM simulations describe the observed 
and M^j distributions reasonably well. 

7.3 Results of the search in the hadronic channel 

The MC simulations of the signal and SM processes were used to find optimal windows 
in M-ij and M^j for the observation of a signal relative to the background. The resulting 
windows were 65.2 < Afjj < 90.8 GeV and 159 < M^J < 188 GeV. 

The closest to M^r is denoted by M^. The distribution of is shown in Fig. El;); 
261 events in the data satisfied the condition 65.2 < < 90.8 GeV. The M^j distri- 
bution after this cut is shown in Fig. Oli). The simulation of SM processes reproduces 
the distributions well. After the requirement 159 < M^^ < 188 GeV, 14 events remained. 
The distributions of and M^j in the data after this cut are shown in Figs. 12^) and f), 
respectively, and are well reproduced by the simulation of SM processes. After these cuts, 
the efficiency for detecting single-top production in the hadronic channel was 24%. This 
efficiency does not include the branching ratio of the top-quark decay in the hadronic 
channel. The observed M^j distribution shows no significant excess at Mtop. 



8 Systematic uncertainties 

The most important sources of systematic uncertainty were: 

• leptonic channel 

— the uncertainty of ±1% on the absolute energy scale of the GAL gave changes of 

in the background and negligible changes in the signal-efficiency estimations; 

— the use of the LEPTO-MEPS model instead of ARIADNE to estimate the NG DIS 
background gave a change of —0.8% in the background estimation; 

— the MG statistical uncertainty on the SM background estimation was ±7.5%; 

• hadronic channel 

— the uncertainty of ±1% on the absolute energy scale of the jets [44] gave changes 
of 1 J°7^% in the background and ^41% in the signal-efficiency estimations; 

— the MG normalisation uncertainty on the SM background estimation was ±7.4%. 
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All these uncertainties in the number of expected background events were added in quadra- 
ture and are shown in Tables ^ and |21 The experimental uncertainties were smaller than 
the theoretical uncertainties and, therefore, were not considered in the derivation of the 
limits for single-top production. 

9 Limit on the FCNC couplings 

As no event was selected in the leptonic channel and no excess over the SM prediction 
was observed in the hadronic channel (see Table |21), limits were set on FCNC couplings of 
the type tqV. The contribution of the charm quark, which has only a small density in the 
proton at high x, was ignored by setting Ktc-y = Vtcz = 0. Only the anomalous couplings 
involving a u quark, Km-^ and v^z^ were considered. 

At HERA, most of the sensitivity to FCNC-induced couplings involving the top quark 
comes from the process ep — etX in which a 7 is exchanged since the large Z mass 
suppresses the contribution due to Z exchange. In a first step, limits on Km-y were, 
therefore, derived assuming Vtuz = and using NLO QCD calculations of the cross section 
for the process ep — > etX (see Section l^)- The results obtained from each channel and 
centre-of-mass energy together with those from the combined analysis presented below are 
summarised in Table |21 Limits from a combination of channels were obtained by using 
a method described in a previous publication [48]. In the derivation of the limits, the 
decrease in the branching ratio B{t ^ Wb) in the presence of FCNC decays was taken 
into account. In comparison to the dependence of the result on the assumed value of Mtop, 
the effects of all other uncertainties are very small. Therefore, limits were evaluated for 
Mtop = 170, 175 and 180 GeV, neglecting the other uncertainties. 

By combining the results from both the leptonic and hadronic channels, an upper limit of 

Ktu^ < 0.174 at 95% CL 

was derived assuming Mtop = 175 GeV. The limit was < 0.158 (0.210) for Mtop = 
170 (180) GeV. The above coupling limit corresponds to a limit on the cross section for 
single-top production of 

a{ep etX, = 318 GeV) < 0.225 pb at 95% CL. 

In a second step, the effects of a non-zero vtuz coupling were taken into account. The 
derivation of the exclusion region in the Htu-y — VtuZ plane was made using LO calculations 
for the process ep etX obtained with the program CompHEP [49], since NLO correc- 
tions to the contribution from Z exchange are not available. Limits in the Kj^i^ — Vtuz 
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plane were derived by using a two-dimensional probability density evaluated assuming a 
Bayesian prior probability distribution flat in ntu-y and Vtuz- 



VtuZ 



D) 



]\iLi{Ni^J\Ktu^,Vtuz) 



Jq dKtu-y Jq dVtuZ Y\iLi{Nl-^J\Ktu-y,VtuZy 



where p^t^tu-y.vtuzlD) is the probability density for the FCNC couplings given the set of 
observed data D and Li{Nl^^\ntu'y^Vtuz) are the partial likelihoods for each channel and 
centre-of-mass energy evaluated as the Poissonian probabilities to observe A^^^j^^ events 
given the expectations of the SM background processes and the signals for single-top 
production. The 95% CL limit was found as the set of points p^K.tu'y.VtuzlD) = po such 
that 



Figure El shows the exclusion region on the Ktu-y — Vtuz plane obtained from this search, 
together with those from CDF [4,50] and L3 [6], which is the most stringent limit from 
LEP2 [5]. It should also be noted that the Lagrangian used in the LEP analyses [5,6,50] 
differs from that in Eq. (1) by a constant multiplicative factor such that nl"^ = \/2 nf^^ 
and = 

v't^z^- In Fig- HI the limits from CDF and L3 are plotted using the 
Lagrangian convention of Eq. (1). The measurements at the Tevatron and LEP have 
similar sensitivities to the tuV and tcV couplings, and their limits were obtained with 
the assumptions = and Vtu-^j = Vtcz- In Fig- HI the published CDF and L3 limits 
are rescaled by \/2 for the purposes of comparison to the present results on Ktu^^ and Vtuz 
which are obtained assuming = Vtcz = 0. The limit-setting procedure was repeated 
assuming Mtop = 170 and 180 GeV; the resulting exclusion regions are also shown in 
Fig.ia 

10 Summary 

Single-top production via flavour- changing neutral current transitions has been searched 
for with the ZEUS detector at HERA in positron-proton and electron-proton collisions at 
centre-of-mass energies of 300 and 318 GeV using an integrated luminosity of 130.1 pb~^. 
No deviation from the Standard Model prediction was found. The results were used to 
constrain single-top production ep etX via the FCNC process. An upper limit on the 
FCNC coupling Km-y of 0.174 at 95% CL was obtained. This limit excludes a substantial 
region in Ktu-y not constrained by previous experiments. 




\D) = 0.95. 
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Leptonic channel 


Positron 
channel 
obs. /expected (H^) 


Muon 
channel 
obs. /expected (W) 




Preselection 


e+p, = 300 GeV {C = 47.9 pb-^) 


4 / 7.3i°i 


/ 4.2t°:^ 


e-p, = 318 GeV (£ = 16.7 ph~^) 


7 / 3.2^?:^ 


1 / 2.11°- 


e+p, = 318 GeV (£ = 65.5 ph~^) 


13 / lO.im 


11 / 5.6i°:l 


Total (£ = 130.1 pb-i) 


24 / 20.6^^:^ (17%) 


12 / 11.9^°:? (16%) 




Final selection [p^-"^ > 25 GeV) 


e+p, ^ = 300 GeV (£ = 47.9 pb-^) 


/ 0.72+°:?37 


/ 0.78l°i° 


e-p, = 318 GeV {C = 16.7 pb-^) 


1 / 0.641°-^ 


1 / 0.45in^ 


e+p, = 318 GeV (£ = 65.5 pb-^) 


1 / 1.54l°:i 


4 / i.53i°:i^ 


Total {C = 130.1 pb-i) 


2 / 2.90+°:! (45%) 


5 / 2.751°;^! (50%) 




Final selection {p^'''^ > 40 GeV) 


e+p, = 300 GeV (£ = 47.9 pb"^) 


/ 0.23+°:°^ 


/ 0.261H1 


e-p, = 318 GeV (£ = 16.7 pb-^) 


0/0.16+°:°^ 


/ 0.08lHi 


e+p, = 318 GeV (£ = 65.5 pb^^) 


/ 0.5At'ofr 


/ 0.61^°;^° 


Total {C = 130.1 pb-i) 


/ OMt'ofo (61%) 


/ 0.95l°:i^ (61%) 



Table 1: Number of events in data and Standard Model background for the 
leptonic channel for different samples after the preselection and final selection cuts. 
The percentage of single-W production included in the expectation is indicated in 
parentheses. The statistical and systematic uncertainties added in quadrature are 
also indicated. 
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Leptonic channel 


Hadronic channel 


300 GeV 


318 GeV 


300 GeV 


318 GeV 










5 


9 




49+0.07 


-'-•^"-'-0.13 


q q+1.3 
0.0_Q 4 




e ■ Br (%) 


6.9 


7.1 


16.6 


16.5 


luminosity (pb^^) 


47.9 


82.2 


45.0 


82.2 


aiim X B{t ^ Wb) (pb) 


0.906 


0.514 


0.998 


0.426 


i^tu-y (per channel) 


0.223 


0.241 


Ciim (pb) (all channels) 


0.225 at = 318 GeV 


Ktu-y (all channels) 


0.174 



Table 2: Number of events in data and Standard Model background for the 
leptonic and hadronic channels for different samples, together with the efficiency 
times branching ratio of the signal and luminosity for each sample. The last four 
rows show the limits on the single-top production cross section via flavour- changing 
neutral current transitions and on the Htu-y coupling assuming M^op = 175 GeV. 
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Figure 1: Single-top production via flavour- changing neutral current transitions 
at HERA. 
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Figure 2: a) $acopj b) p^'^, and c) for those events with an identified 
positron candidate, d) p^^ , e) p^'^ and f) pi^ for those events with an identified 
muon candidate. The dots are the data, the solid histogram is the Standard Model 
MC simulation and the shaded histogram represents the signal with Mtop = 175 GeV 
normalised to the limit presented in Section\^ The final bins in d) and f), marked 
"overflow" , contain all events above the lower boundaries of these bins. The dis- 
tributions are for the selected events according to the criteria of Section \6. 1\ The 
Standard Model MC distributions have been normalised to the luminosity of the 
data. 
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ZEUS 





Figure 3: a) M-^j (for all pairs of jets), b) M^j and c) M^^ distributions for 
the sample of events selected in the hadronic channel; d) M^^ distribution for those 
events with 65.2 < My^ < 90.8 GeV; e) My^ distribution for those events with 159 < 
M^j < 188 GeV; f) M^J distribution for those events with 65.2 < < 90.8 GeV. 
The dashed lines in c) and e) represent the cut at 65.2 < My^ < 90.8 GeV. The 
distributions are for three-jet events with E^'^ > 40, E^^*^ > 25, E^^'^ > U GeV 
and —1 < rf^^ < 2.5. Other details are given in the caption to Fig. 
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Figure 4: Exclusion regions at 95% CL in the Htu-y — v^z plane for three values 
of Mtop (170, 175 and 180 GeV) assuming ntc-y = Vtcz = 0. The CDF and L3 
exclusion limits for Mtop = 175 GeV are also shown. 
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Limits using an improved signal simulation 



In this addendum, a new simulation of single-top production, improved with respect to 
that described in Section 4, has been employed. A full simulation of the flavour-changing 
neutral current process has been performed, considering the anomalous tuZ and tw^ 
couplings {vtuz^ f^tu-y) both in the top production and decay. This allows the exclusion 
limits in the Ktwy — Vtuz plane to be improved at large values of Vtuz- 

The exclusion limits obtained in Section 9 used the HEXF generator [1], which considers 
only 7-mediated top production to evaluate the selection efficiency for the signal. The 
effect of the Z-mediated process was taken into account by the conservative assumption 
that the efficiency was zero in the hadronic channel and identical to that for the 7 case in 
the leptonic channel. However, the final states from top production mediated by 7 or Z 
exchange are quite different. For the former, the polar- angle distribution of the scattered 
electron^ is peaked towards the electron beam direction and often escapes detection, while 
for the latter, the scattered electron is within the detector acceptance in most cases. As 
a consequence, the efficiency for the Z-exchange process is larger than that for the 7- 
exchange process in the leptonic decay channel, and lower in the hadronic decay channel, 
where an explicit cut against final-state electrons was applied to suppress the neutral 
current DIS background. 

The improved simulation uses COMPHEP [2] for the hard process eu —>■ et and 
PYTHIA6 [3] for proper treatment of the proton remnant, initial- and final-state QED and 
QCD radiation, hadronisation and particle decays. Monte Carlo samples were generated 
for the top decay channels t — > bW and t uZ. For 7 exchange, HEXF and COMPHEP 
give consistent results. Final results have been obtained using HEXF for 7 exchange and 
COMPHEP for Z exchange. The interference between the two processes gives a contri- 
bution to the cross section of less than 1% in the whole range of the anomalous couplings 
considered in the analysis and was neglected. 

For the process mediated by a photon, the selection efficiency in the leptonic channel 
is 7.1% (3.0%) for t — > hW (t — ^ uZ), while for the hadronic channel it is 16.5% for 
t bW and slightly higher for t — » uZ. The new simulation allows the evaluation of 
the efficiency for the Z-exchange process. In the leptonic channel, it is 11.5% (10.6%) for 
t — > bW {t —>■ uZ), while for the hadronic channel it is 4.0% for both decay channels. The 
reported efficiencies are for ^/s = 318 GeV and include the branching ratio of the W or 
Z in the considered decay channels; for ^/s = 300 GeV, the efficiencies are similar. The 
efficiency of the selection criteria, described in Section 6.7, is small for the decay channel 
t —>■ wf, and was conservatively assumed to be zero. 

^ Here and in the following 'electron' is used indifferently for both electron and positron. 
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The derivation of the exclusion region in the Ktu-y — Vtuz plane was made using LO cal- 
culations for the process ep — > etX obtained with the program COMPHEP. Corrections 
due to QED initial-state radiation were taken into account using the Weizsacker- Williams 
approximation [4,5]. Figure Q shows the 95% confidence level exclusion limits in the 
Ktu^ — Vtuz plane obtained from this search, together with those from CDF [6], L3 [7] and 
HI [8]. The HI result constrains only the ntu-y coupling; it was obtained including next- 
to- leading-order QCD corrections and assuming Vtuz = 0. The limits were obtained using 
the procedure described in Section 9. The improved simulation of single-top production 
considerably improves the limits for large values of the anomalous vector coupling Vtuz- 
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Figure 1: Exclusion regions at 95% CL in the Htuj — Vtuz plane for three values 
o/Mtop (170, 175 and 180 GeV) assuming Ktc-^ = vtcz = 0. The CDF, L3 and HI 
exclusion limits are also shown. 
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